The conserved Arg241‐Glu439 salt bridge determines flexibility of the inositol 1,4,5‐trisphosphate receptor binding core in the ligand‐free state

Inositol 1,4,5‐trisphosphate receptor (InsP₃R) is an intracellular Ca²⁺‐release channel activated by binding of inositol 1,4,5‐trisphosphate (InsP₃) to the InsP₃ binding core (IBC). Structural change in the IBC upon InsP₃ binding is the key process in channel pore opening. In this study, we performed molecular dynamics (MD) simulations of the InsP₃‐free form of the IBC, starting with removal of InsP₃ from the InsP₃‐bound crystal structure, and obtained the structural ensemble of the InsP₃‐free form of the IBC. The simulation revealed that the two domains of the IBC largely fluctuate around the average structure with the hinge angle opened 17° more than in the InsP₃‐bound form, and the twist angle rotated by 45°, forming interdomain contacts that are different from those in the bound form. The InsP₃ binding loop was disordered. The InsP₃‐free form thus obtained was reproduced four times in simulations started from a fully extended configuration of the two domains. Simulations beginning with the fully extended form indicated that formation of a salt bridge between Arg241 and Glu439 is crucial for stabilizing the closed form of the two domains. Mutation of Arg241 to Gln prevented formation of the compact structure by the two domains, but the fully flexible domain arrangement was maintained. Thus, the Arg241‐Glu439 salt bridge determines the flexibility of the InsP₃‐free form of the IBC.Proteins 2013; 81:1699–1708. © 2013 Wiley Periodicals, Inc.     

Deriving how far structural information is transmitted through parallel homodimeric coiled‐coils: A correlation analysis of helical staggers

How local conformation is affected by local sequence is fairly well understood for alpha‐helical coiled‐coils, but less is known about how local conformation is influenced by distant features. Here, I describe an approach to detect such an effect, based on computing correlation coefficients of local out‐of‐register alignments, or so‐called “staggers” between the helices, as a function of the axial distance between the staggers. This approach requires parallel homodimers, in which each stagger can occur with two “signs,” where either one helix or the other is shifted towards the N terminus. The signs of such staggers separated by up to 12 residues are strongly correlated, indicating that the conformations of the ends of coiled‐coils are commonly influenced by attached structures. Thus, the structures of coiled‐coil residues aberrantly attached to alternative proteins, such as those resulting from leukemogenic chromosomal rearrangements, may be distinguishable from those in normal tissues, and in turn serve as targets of selective drug design. The signs of helical staggers separated by between 13 and 30 residues are moderately yet significantly correlated, indicating that some of the coiled‐coils transmit this conformational feature axially for at least 45 Å. A positive, albeit noisy, correlation also exists among tropomyosin coiled‐coils for signed staggers separated by the 40‐residue actin repeat distance, consistent with the semi‐flexible tropomyosin filament binding F‐actin and regulating skeletal muscle contraction in a partially cooperative manner. Communication of the signs of axial staggers is explained in part by minimization of main‐chain hydrogen bond deformations. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Sequence and hydropathy profile analysis of two classes of secondary transporters

A structural class in the MemGen classification of membrane proteins is a set of evolutionary related proteins sharing a similar global fold. A structural class contains both closely related pairs of proteins for which homology is clear from sequence comparison and very distantly related pairs, for which it is not possible to establish homology based on sequence similarity alone. In the latter case the evolutionary link is based on hydropathy profile analysis. Here, we use these evolutionary related sets of proteins to analyze the relationship between E-values in BLAST searches, sequence similarities in multiple sequence alignments and structural similarities in hydropathy profile analyses. Two structural classes of secondary transporters termed ST[3], which includes the Ion Transporter (IT) superfamily and ST[4], which includes the DAACS family (TC# 2.A.23) were extracted from the NCBI protein database. ST[3] contains 2051 unique sequences distributed over 32 families and 59 subfamilies. ST[4] is a smaller class containing 399 unique sequences distributed over 2 families and 7 subfamilies. One subfamily in ST[4] contains a new class of binding protein dependent secondary transporters. Comparison of the averaged hydropathy profiles of the subfamilies in ST[3] and ST[4] revealed that the two classes represent different folds. Divergence of the sequences in ST[4] is much smaller than observed in ST[3], suggesting different constraints on the proteins during evolution. Analysis of the correlation between the evolutionary relationship of pairs of proteins in a class and the BLAST E-value revealed that: (i) the BLAST algorithm is unable to pick up the majority of the links between proteins in structural class ST[3], (ii) ‘low complexity filtering’ and ‘composition based statistics’ improve the specificity, but strongly reduce the sensitivity of BLAST searches for distantly related proteins, indicating that these filters are too stringent for the proteins analyzed, and (iii) the E-value cut-off, which may be used to evaluate evolutionary significance of a hit in a BLAST search is very different for the two structural classes of membrane proteins.     

Relative roles of dispersal dynamics and competition in determining the isotopic niche breadth of a wetland fish

1. The niche variation hypothesis predicts that among‐individual variation in niche use will increase in the presence of intraspecific competition and decrease in the presence of interspecific competition. We sought to determine whether the local isotopic niche breadth of fish inhabiting a wetland was best explained by competition for resources and the niche variation hypothesis, by dispersal of individuals from locations with different prey resources or by a combination of the two. We analysed stable isotopes of carbon and nitrogen as indices of feeding niche and compared metrics of within‐site spread to characterise site‐level isotopic niche breadth. We then evaluated the explanatory power of competing models of the direct and indirect effects of several environmental variables spanning gradients of disturbance, competition strength and food availability on among‐individual variation of the eastern mosquitofish (Gambusia holbrooki). 2. The Dispersal model posits that only the direct effect of disturbance (i.e. changes in water level known to induce fish movement) influences among‐individual variation in isotopic niche. The Partitioning model allows for only direct effects of local food availability on among‐individual variation. The Combined model allows for both hypotheses by including the direct effects of disturbance and food availability. 3. A linear regression of the Combined model described more variance than models limited to the variables of either the Dispersal or Partitioning models. Of the independent variables considered, the food availability variable (per cent edible periphyton) explained the most variation in isotopic niche breadth, followed closely by the disturbance variable (days since last drying event). 4. Structural equation modelling provided further evidence that the Combined model was best supported by the data, with the Partitioning and the Dispersal models only modestly less informative. Again, the per cent edible periphyton was the variable with the largest direct effect on niche variability, with other food availability variables and the disturbance variable only slightly less important. Indirect effects of heterospecific and conspecific competitor densities were also important, through their effects on prey density. 5. Our results support the Combined hypotheses, although partitioning mechanisms appear to explain the most diet variation among individuals in the eastern mosquitofish. The results also support some predictions of the niche variation hypothesis, although both conspecific and interspecific competition appeared to increase isotopic niche breadth in contrast to predictions that interspecific competition would decrease it. We think this resulted from high diet overlap of co‐occurring species, most of which consume similar macroinvertebrates.     

A 5q12.1–5q12.3 microdeletion in a case with a balanced exceptional complex chromosomal rearrangement

Complex chromosomal rearrangements are very rare chromosomal abnormalities. Individuals with a complex chromosomal rearrangement can be phenotypically normal or display a clinical abnormality. It is believed that these abnormalities are due to either microdeletions or microduplications at the translocation breakpoints or as a result of disruption of the genes located in the breakpoints. In this study we describe a 2-year-old child with mental retardation and developmental delay in whom a de novo apparently balanced exceptional complex chromosomal rearrangement was found through conventional cytogenetic analysis. Using both cytogenetic and FISH analysis, the patient's karyotype was found to be: 46,XY,der(5)t(5;7)(p15.1;7q34),t(5;8)(q13.1;8q24.1)dn. A large, clinically significant deletion which encompassed 887.69 kb was detected at the 5q12.1–5q12.3 (chr5:62.886.523–63.774.210) genomic region using array-CGH. This deleted region includes the HTR1A and RNF180 genes. This is the first report of an individual with an apparently balanced complex chromosomal rearrangement in conjunction with a microdeletion at 5q12.1–5q12.3 in which there are both mental-motor retardation and dysmorphia.     

Characterization of the first adult de novo case of 46,X,der(Y)t(X;Y)(p22.3;q11.2)

Herein, we describe a case of an infertile man detected in postnatal diagnosis with FISH characterization and array-CGH used for genome-wide screening which allowed the identification of a complex rearrangement involving sex chromosomes, apparently without severe phenotypic consequences. The deletion detected in our patient has been compared with previously reported cases leading us to propose a hypothetical diagnostic algorithm that would be useful in similar clinical situations, with imperative multi disciplinary approach integrated with genetic counseling. Our patient, uniquely of reproductive age, is one of six reported cases of duplication of Xp22.3 (~ 8.4 Mb) segment and contemporary deletion of Yq (~ 42.9 Mb) with final karyotype as follows:

46,X,der(Y),t(X;Y)(Ypter → Yq11.221::Xp22.33 → Xpter).ish der(Y) (Yptel+,Ycen+,RP11-529I21+,RP11-506M9-Yqtel −,Xptel +). arrXp22.33p22.31(702–8,395,963, 8,408,289x1), Yq11.221q12 (14,569,317x1, 14,587,321–57,440,839x0)     

Comparing methods for metabolic network analysis and an application to metabolic engineering

Bioinformatics tools have facilitated the reconstruction and analysis of cellular metabolism of various organisms based on information encoded in their genomes. Characterization of cellular metabolism is useful to understand the phenotypic capabilities of these organisms. It has been done quantitatively through the analysis of pathway operations. There are several in silico approaches for analyzing metabolic networks, including structural and stoichiometric analysis, metabolic flux analysis, metabolic control analysis, and several kinetic modeling based analyses. They can serve as a virtual laboratory to give insights into basic principles of cellular functions. This article summarizes the progress and advances in software and algorithm development for metabolic network analysis, along with their applications relevant to cellular physiology, and metabolic engineering with an emphasis on microbial strain optimization. Moreover, it provides a detailed comparative analysis of existing approaches under different categories.     

Structural analysis of the envelope gp120 V3 loop for some HIV-1 variants circulating in the countries of Eastern Europe

The V3 loop on gp120 from human immunodeficiency virus type 1 (HIV-1) is a focus of many research groups involved in anti-AIDS drug development because this region of the protein is a principal target for neutralizing antibodies and a major determinant for cell tropism and syncytium formation. In this study, the nucleotide sequences of the env gene region coding the V3 loop were determined by DNA sequencing methods for four novel HIV-1 strains that circulate in the countries of Eastern Europe, such as Russia, Belarus, Ukraine, etc. Based on the empirical data obtained, the 3D structures of the V3 loops associated with these viral modifications were generated by computer modeling and then compared to discover similarities in the spatial arrangement of this functionally important site of gp120. Despite the HIV-1 genetic variety, several regions of the V3 loop that contain residues critical for cell tropism were shown to be structurally invariant, which may explain its exceptional role in a co-receptor usage. These data together with those on the biological activity of the V3 individual residues clearly show that these conserved structural motifs of gp120 represent potential HIV-1 weak points most suitable for therapeutic intervention.     

Did the prion protein become vulnerable to misfolding after an evolutionary divide and conquer event?

Despite high sequence identity among mammalian prion proteins (PrPs), mammals have varying rates of susceptibility to prion disease resulting in a so-called species barrier. The species barrier follows no clear pattern, with closely related species or similar sequences being no more likely to infect each other, and remains an unresolved enigma. Variation of the conformationally flexible regions may alter the thermodynamics of the conformational change, commonly referred to as the conformational conversion, which occurs in the pathogenic process of the mammalian prion protein. A conformational ensemble scenario is supported by the species barrier in prion disease and evidence that there are strains of pathogenic prion with different conformations within species. To study how conformational flexibility has evolved in the prion protein, an investigation was undertaken on the evolutionary dynamics of structurally disordered regions in the mammalian prion protein, non-mammalian prion protein that is not vulnerable to prion disease, and remote homologs Doppel and Shadoo. Structural disorder prediction analyzed in an evolutionary context revealed that the occurrence of increased or altered conformational flexibility in mammalian PrPs coincides with key events among PrP, Doppel, and Shadoo. Comparatively rapid evolutionary dynamics of conformational flexibility in the prion protein suggest that the species barrier is not a static phenomenon. A small number of amino acid substitutions can repopulate the conformational ensemble and have a disproportionately large effect on pathogenesis.     

Molecular characterization influencing metal resistance in the Cupriavidus/Ralstonia genomes

Our environment is stressed with a load of heavy and toxic metals. Microbes, abundant in our environment, are found to adapt well to this metal-stressed condition. A comparative study among five Cupriavidus/Ralstonia genomes can offer a better perception of their evolutionary mechanisms to adapt to these conditions. We have studied codon usage among 1051 genes common to all these organisms and identified 15 optimal codons frequently used in highly expressed genes present within 1051 genes. We found the core genes of Cupriavidus metallidurans CH34 have a different optimal codon choice for arginine, glycine and alanine in comparison with the other four bacteria. We also found that the synonymous codon usage bias within these 1051 core genes is highly correlated with their gene expression. This supports that translational selection drives synonymous codon usage in the core genes of these genomes. Synonymous codon usage is highly conserved in the core genes of these five genomes. The only exception among them is C. metallidurans CH34. This genomewide shift in synonymous codon choice in C. metallidurans CH34 may have taken place due to the insertion of new genes in its genomes facilitating them to survive in heavy metal containing environment and the co-evolution of the other genes in its genome to achieve a balance in gene expression. Structural studies indicated the presence of a longer N-terminal region containing a copper-binding domain in the cupC proteins of C. metallidurans CH3 that helps it to attain higher binding efficacy with copper in comparison with its orthologs.